Scattering of atoms on a Bose— Einstein condensate 
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We study the scattering properties of a Bose-Einstein condensate held in a finite depth well when 
the incoming particles are identical to the ones in the condensate. We calculate phase shifts and 
corresponding transmission and reflection coefficients, and we show that the transmission times can 
be negative, i.e., the atomic wavepacket seemingly leaves the condensate before it arrives. 

PACS numbers: 3.75.Fi,34.50-s 
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I. INTRODUCTION 

In this paper we present an analysis of the scattering of 
atoms on a Bose-Einstein condensate. A repulsive inter- 
action between atoms is both responsible for the shape 
of the condensate and for the interaction between the in- 
cident atom and the condensed atoms. Scattering effects 
can be significant when two condensates merge, and four- 
wave mixing has been observed in experiments We 
deal here with the limit of a single or just a few atoms, 
incident in a well defined momentum state on a conden- 
sate, a situation that has been realized experimentally to 
seed atom lasers Q and 'atomic parametric amplifiers' 
|5|, |^| with a weak atomic beam. Our interest, however, 
is in the interaction between the small atomic compo- 
nent and the condensate, and in particular in the trans- 
mission and reflection properties of the condensate as a 
beam splitter for atoms of the same kind as the mirror 
itself. Since the atoms of the scatterer are indistinguish- 
able from the scattered particles, exchange effects play 
an important role, and similar scattering studies have 
indeed been proposed as a means for investigatingsuper- 
fluidity in strongly interacting He-4 system ||, f|, and 
also in weakly interacting systems . 

Dilute condensates can be trapped in a large number 
of trapping arrangements, and we shall assume a trap- 
ping potential of finite width allowing asymptotically free 
atoms to be directed towards the trapped condensate and 
interact with it in a well defined region of space. Such 
potentials can be created optically or possibly in the re- 
gion above a current carrying chip, and our calculations, 
which will be carried out in a one dimensional geometry, 
may correspond to either a condensate slab hit by atoms 
with well defined transverse momenta or to a trapped 
condensate in a wave guide with a local longitudinal min- 
imum. 

In Sec. II we introduce the Bogoliubov treatment of 
excitations of a Bose-Einstein condensate. In Sec. Ill we 



explain how an analysis of the fundamental excitations of 
the condensate provides the scattering information. Nu- 
merical results are compared with an approximate an- 
alytical model. In Sec. IV, we show examples of time 



dependent wavepacket dynamics, and we identify a mo- 
mentum regime in which transmission occurs with nega- 
tive time delays, i.e., the wave packet emerges from the 
condensate before it arrives. Sec. concludes the paper. 



II. THE BOGOLIUBOV APPROACH 

In the Bogoliubov approach, the number of particles 
in the condensate mode is assumed to be high and rather 
well defined, and the problem of small thermal or ex- 
ternally induced excitations can be treated to lowest or- 
der to yield a picture of non-interacting quasi-particles. 
These are the normal modes of the system and they 
are mixtures of particle- and hole-like excitations. In 
the case of scattering, the incoming particle will behave 
as a quasiparticle while moving through the condensate, 
whereas asymptotically the excitation is particle-like as 
hole-excitations are limited to the region occupied by the 
condensate. 

Let us briefly review the Bogoliubov approach. First, 
the field operator is written in a form explicitly empha- 
sizing the condensate mode 



W(x) = VN(j} (x) + 5^(x). 



(1) 



Here 4>o(x) is to be regarded as a c-number field de- 
scribing the condensate while 6^(x) describes the cor- 
rection. We insert this form in the full Hamiltonian (in 
the contact-interaction approximation), 



dx \ &(x)h(x)i>(x) 



H 



+ ^* t (x)* t (x)*(a;)*(x)|, (2) 

where h(x) = — (h 2 /2m)d 2 + V(x) is the single-particle 
Hamiltonian and giD quantifies the ID interaction 
strength. Keeping only interaction terms of order N 
or higher, we arrive at a quadratic form in S^(x) and 
S^f'(x). The terms linear in ^(x) and ^(x) are identi- 
cally zero if (f>o(x) solves the Gross-Pitaevskii equation 
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[h{x) + gN\Mx)\ 2 ] <t>o{x) = M0o(»- 



(3) 



2 



To diagonalize the remaining part, one envokes a Bogoli- 
ubov transformation and writes 



S9(x) = [ku k (x) - b\v* k {x) 



(4) 



where Uk(x) and v k (x) solve the Bogoliubov-de Gcnncs 
equations: 

[L(x) - hujk]uk(x) = gn a (x)v k {x) (5) 
[L{x) + hujk]vk(x) = gn a (x)u k (x) (6) 

and arc normalized to 



\u* k (x)uk'{x) - v%(x)vk'{x)]dx = 5 k 



k 1 ■ 



(7) 



The operator L(x) is defined by L(x) = h(x) + 2gno(x) 
jj,. The Hamiltonian finally takes the simple form 



Ha 



E(N) 



OO 

E 

k=l 



Tiu} k b\b k , 



(8) 



i.e, it describes non-interacting quasi-particles created 
by b\ and destroyed by b k - Note how the Bogoliubov 
transform mixes creation and destruction operators. This 
means, that excitations have both a particle-like charac- 
ter (described by Uk) and a hole-like character (describe 
by v k ). 

The Bogoliubov analysis is routinely applied to the 
case of a homogeneous situation (flat potential) and the 
case of an infinite trapping potential In the first 
case, one obtains solutions with well-defined momenta: 
phonons in the long wave-length limit, free particles for 
high momenta. The spectrum is continuous and gap-less. 
In the second case, the excitations are again collective in 
the low energy regime, approaching single-particle trap 
states for high energies ||. The spectrum is discrete. 
Here we will concentrate on a third situation where the 
trapping potential has finite width and depth. In this 
scenario, a finite number of trapped excitations exists 
and above these a continuum of scattering states. The 
corresponding spherically symmetric problem in 3D has 
been treated by Wynveen et al. in . 

Away from the trap and thus from the condensate, the 
scattering states are particle-like (vk(x) = 0) and the 
usual asymptotic analysis of scattering applies: We can 
define incoming, reflected, and transmitted fluxes and we 
can set up an initial wave-packet and make it progate to- 
wards the condensate region. The actual tranmission and 
reflection coefficients depend on what happens while the 
incoming particle passes through the condensate and this 
propagation is partly phonon like, i.e., with the emer- 
gence of a hole-like component. 

In ID, the Bogoliubov-De Gennes equations are com- 
putationally quite manageable and the choice of method 
for their solution is not crucial. Nevertheless, for com- 
pleteness we here briefly summarize our method. First 
we solve the Gross-Pitaevskii equation (||) with the cho- 
sen potential. We do this by a simple steepest descent 
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FIG. 1: Figure showing the potential well and the shape of 
the condenste. The range of the potential is [—a, a] outside 
which V(x) is identically 0. V(x) changes smoothly to — Vo = 
— 250H 2 /ma 2 in edge-zones of width 0.3a. The ID interaction 
strength is taken to be gmN — 400ft 2 /ma. The solution of 
the GPE (|) then gives /x = - 14. 4h 2 /ma 2 and a Bogoliubov 



speed of sound c = sj gno(0)/m 
region of the well. 



15. Ah /ma in the inner 



method, i.e., by propagating the corresponding time- 
dependent equation in imaginary time. In particular, we 
use a split-step fast-Fourier-transform algoritm |T7J. In 
Fig. [j], we show the particular potential we will be using, 
and we show the corresponding condensate wavefunction. 

When 4>o(x) and /i have been found, Eqs. (|5|) and (|^) 
are solved by first defining 



(9) 



fk(x) = V g i U k( X ) + v k{x)) 



hk{x) = \j - [u k (x) - v k (x) 



in terms of which Eqs. (j5|) and (||) become 

[L(x) - gn (x)]f k (x) = huj k h k {x) 
[L(x) + gn (x)}h k (x) = huj k fk{x) 



(10) 
(11) 



We can then avoid diagonalizing a two-component prob- 
lem by simply applying the operator [L(x) + gno(x)] to 
both sides of Eq. ( |l0|) to find the necessary condition 

[L(x)+gn (x)][L(x)-gno(x)]f k (x) = [hco k ] 2 f k (x). (12) 

When this one-component eigenvalue problem is solved, 
Uk(x) and v k {x) can be found by applying first Eq. jlO| ) 
and then Eq. (ft). 
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III. RESULTS 



A. Phaseshifts 



We assume that V(x) and thus no{x) have even spa- 
tial symmetry. We can then demand that solutions of 
Eqs. <JsJ) and (||) have definite parity. In the discrete 
part of the spectrum, even and odd solutions alternate, 
while in the continuum, each energy supports solutions 
of both parities. In the asymptotic region away from the 
potential and the condensate, even and odd scattering 
solutions can then be written 



4 e) 
» 



(x) — > cos(fca; =F 5 e (fc)) 
(x) —>■ sin(kx =F <5 D (fc)) 



for x — > ±00. 



(13) 



Note that we identify the label k with the asymptotic 
wavenumber for the scattering solutions. The phase- 
shifts, S e (k) and <5 c (fc), contain all the information about 
the scattering relevant to the asymptotic region. They 
can be conveniently extracted from our numerical solu- 
tions for Uk- The calculations are actually done on a 
space-interval [— L, L] with periodic boundary conditions 
(finite lower momentum cutoff). This means that for a 
given L we only find the subset of the continuum solu- 
tions with wavenumbers fulfilling 



-kL + 6(k) = kL-6(k) mod 2tt 



or, equivalcntly, 



S(k) = kL mod tt. 



(14) 



(15) 



After the diagonalization this discrete set of k values 
can easily be determined from the eigenvalues as tuvk = 
h 2 k 2 /2m — fx. Each fc value gives us one point on 5(k) 
via Eq. (|l5|). This is enough to determine S(k) if it is 
slowly varying on the 1/L scale. If not, we just need to 
change L slightly and repeat the calculation to obtain an 
additional set of points. 

In Fig. |^ we show a typical example of S e (k) and S (k) 
with parameters as in Fig. |IJ There is clearly some reso- 
nant behavior with out-of-phase oscillations of 5 e (k) and 
6 (k) around a slowly varying average. This behaviour 
is well known from ordinary, single particle scattering on 
a well/barrier, and in the following subsection we shall 
present an analytical model which yields the same gross 
features. 



B. Square well model, Thomas— Fermi 
approximation 




FIG. 2: Phaseshifts of even and odd solutions as a function 
of incoming momentum. Parameters as in Fig. [l]. 



where a aw is an appropriately defined effective width. In 
the Thomas-Fermi approximation the condensate wave 
function is constant in the trap and zero outside 



(j) (x) 



, |x| > a s 



(17) 



This 4>o(x) is naturally not a solution to the Gross- 
Pitaevskii equation as the kinetic energy term will 
smoothen the step in density even when the external po- 
tential is discontinuous. Formally, this introduces linear 
terms in Eq. (|sj> , i.e., the Bogoliubov vacuum is not a 
steady-state for the system. We are, however, only in- 
terested in the scattering behaviour at positive energies 
and it is reasonable to assume that some insight can be 
gained by finding solutions to Eqs. (||) and (g[) with the 
simplifying assumptions expressed by ( |l6| ) and (p^). 

It is amusing to note, that we are now dealing with 
the Bogoliubov-de Gennes analog of the undergraduate 
textbook problem of ID scattering on a square well. The 
full solution is found by matching the analytical wave- 
functions in the regions, x € [— 00, — Osw]j [— <JSW) ssw], 
and [osW) 00 ]- To the left and to the right of the well, 
we demand v(x) — and let u(x) — cos(kx ± 5 e (k)) 
(u(x) = sin(fcx ± (5 Q (fc))) to find even (odd) solutions. 
Inside the well, we are in a region of constant poten- 
tial (V(x) = —Vo) and constant condensate density 
{N\Mx)\ 2 = n o) so Eqs. (§) and (§) read 



Let us consider a condensate trapped in a square well 
potential 



, -Vo , \x\ < a sw 
V( - X) -^ , \x\>a sw 



(16) 



~ — dl - V + 2gn - fi-hu; 



2m 



d 2 — Vo + 2gn a - [i + Tiuj 



u(x) = gn v(x) (18) 
v(x) = gnou(x). (19) 
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To be able to match the boundary conditions, we need 
four linearly independent solutions and the ansatz u{x) = 
Ue lXx , v(x) = Ve lXx leads to A <E {±A,±A} where 

h\(k) = 

J 2m J^l [gn ] 2 + [E - fi] 2 - 2gn + V + ^ (20) 



and 



hX(k) 



ul2m [ \/[gn Q } 2 + [E - ^] 2 + 2gn -Vq— fi) (21) 



where E = (hk) 2 /2m is the incoming kinetic energy. 
Note that A is imaginary and it is normally disregarded 
in homogeneous condensates. Here, it is needed as we 
must match the boundary conditions that also v and v' 
are continuous. 

The matching of u and v at x — a sw leads to equations 
for the phaseshifts. They read 



tan (ka sw — S e (k)) = 
tan (ka sw — S (k)) = 



A(fc) 
k 
k 



tan (A(fc)a s 



(22) 



tan (A(fe)a sw ) . (23) 



These equations are the same as for single pa rticle square 
well scattering except that the usual hX = \j2m(E — V) 
is replaced by TlX from Eq. (^o|). 

In Fig. ^ we show results obtained for parameters like 
in Fig. Ec grig is taken as the central density in the smooth 
trap and a sw is chosen to accomodate the same total 
number of particles. A comparison with Fig. reveals 
both similarities and differences: The smooth behaviour 
of the average of the two curves is well reproduced as 
well as the (quasi-)period of the oscillatory behaviour. 
However, the positions of curve crossings and the maxima 
of the phaseshift differences, especially at low k's, are not 
well reproduced. Also, at high momenta the sharp-edge 
approximation leads to stronger oscillations in S e (k) and 
S (k) than for the smooth well. 



7 
6 
5 
4 
3 
2 
1 

-1 







r\ A 




■-■S 




/ / 


/ 








/ 1 
1/ 








f 1 










•ft 
/ / 










// 

























10 15 20 

ka 



25 



30 



FIG. 3: Even and odd phaseshifts for a simple square 
well/square condensate model. The parameters are like in 
Fig. and the curves should be compared to the full numerical 
calculations presented in Fig. ^. There is clearly a qualitative 
agreement: The smooth behavior of the average of the two 
curves is well reproduced and the periode of the oscillations 
around this average is also of the correct order of magnitude. 
However, the oscillations extend to too high values of k. 



and even parity eigenstates (|l3|), one finds 

= j e -*(*.(*)+f.W) sin(<5 (fc) - 5 e (k)) 
T(jfe) = i(e- 2i4 'W+e- 2M -W) 



= e - J (^( fe )+ 5 =W) cos(,5 (fc) - 6 e (k)). 



(25) 



(26) 



We see that 100% transmission takes place at fc's where 
S e (k) — 5 (k) while 100% reflection requires 5 e (k) = 
5 (k) ± 7r/2. In Fig. || we plot R and T for the same 
situation as considered above. The curve crossings in 
Fig. |^ are now translated into transmission windows. At 
very low momenta, we get 100% reflection and at very 
high momenta we naturally get 100% transmission. 



Transmission coefficent 



IV. TIME DEPENDENT SCATTERING 



The typical scattering situation with an incoming, a re- 
flected, and a transmitted wave specifies the asymptotic 
form of the wave function 



, , J e lkx + R(k)e lkx for x — > -co 
Uk{x) | T ^y kx for x ^ OQ 



(24) 



This defines the reflection- and transmission-coefficients 
R(k) and T(k). Making the change of basis from the odd 



A. Wavepacket dynamics 

With the complete set of scattering states it is possible 
to follow scattering of wavepackets in time. To obtain the 
initial state we add a number of particles to the Bogoli- 
ubov vacuum: 



\t = 0) = 



1 



U^)S^(x) 



vac) 



(27) 
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FIG. 4: Transmission and reflection coefficients for the situ- 
ation considered in Figs, [j] and ^. A series of perfect trans- 
missions are seen, corresponding to the crossings of the S e (k) 
and So(k) curves in Fig. H. 



where 4> s (x) is the desired wavepacket mode. The oper- 
ator term in this equation creates an n s -particle state, 
but to benefit from the simple form of Eq. (|^) we should 
think of it as creating n s quasiparticles which at t = 
just happen to be localized well away from the conden- 
sate. The quasiparticles arc created in a superposition of 
energy eigenstates and the coefficients in this superposi- 
tion are found as 



Cfe 



K(x)0 s (x) + vl(x)(/)*{x)} dx. 



(28) 



As <f> s is located well away from the condensate region, 
there is in fact no contribution from the v part of this 
integral. 

The time evolution is entirely given by the relation 
bk(t) — e~ ZUJkt bk in the Heisenberg picture, and we find 
the non-condensate mode part of the total density to be 



(^(a:)** (a:)) = n s \U(x,t)\ 2 +n s \V(x,t)\ 2 +Y^ I 

k 

where 



U(x,t)=Y,c k e-^ t u k (x) 

k 

V{x,t)=Y^c k e-^ t v k {x). 



(29) 

(30) 
(31) 



The last term in Eq. (|29| ) is the quantum depletion, which 
is always present due to the interactions in the conden- 
sate, while the first two terms are consequences of the 
scattering process. 



Depending on the spread of fc-values in the wavepacket, 
the asymptotic behaviour of the scattering will be more 
or less simply described by reading off the reflection and 
transmission coefficients R and T at the average momen- 
tum fco- In fact, looking at the first form of these coef- 
ficients in Eqs. ( p5|) and ((2^) we are reminded that the 
reflected and the transmitted wave can be seen as super- 
positions of an even part and an odd part. To each of 
these can be ascribed a time-delay, At e / in the arrival 
of the original wave packet at certain point in space with 
respect to the free propagation. It is easy to show that 
these delays are given by 



At n 



2 dS e (k) 



v g dk 
2 dS Q {k) 



Ok 



k=ko 



k—kn 



(32) 
(33) 



where v g — hko jra is the (group) velocity of the incoming 
wavepacket. If these delays are sufficiently different the 
reflected and transmitted wavepackets are expected to be 
doublepeaked. In fact, a closer look at Eqs. ( p5| ) and ( |3"3| ) 
reveals a difference in sign: the reflected wavepacket can 
easily be double-peaked as it is the difference between 
the even and odd contribution, while the transmitted 
wavepacket is the sum and therefore require displace- 
ments as large as the wavepacket width for a visible effect. 

In Fig. |5j we show time series of snapshots from 
wavepacket simulations. The double peak phenomenon 
mentioned above is visible in the k^a = 11.20 timeserie. 
This ko corresponds to a crossing of 6 e (k) and 5 (k) and 
thus to both high tranmission and a marked difference in 
At e and At (cf. Fig. |). 

B. Transmission times 

The time-delays of Eqs. ([52|) and ( |33|) can also be 
translated into an effective transmission time, the time 
spent traversing the condensate region. The well has a 
width 2a so the time spent inside the well can heuristicly 
be defined as r = At + 2a /v g . In Fig. || we plot r e (fc) 
and r (fc). The two curves agree for high k, while for 
ka less than ~ 15, i.e., for v g less than the Bogoliubov 
speed of sound in the homogeneous part of the conden- 
sate, c = y/gno/m, they show alternating peaks. Such 
peaks are signatures of resonances where an even or an 
odd number of oscillations fit inside the condensate. 

At low fc-values, we observe that r e becomes nega- 
tive over a rather wide range. For wavepackets with mo- 
mentum components mainly in this range, a peak in the 
transmitted wavepacket can appear before the peak of the 
incident wavepacket has reached the condensate. This is 
confirmed by wavepacket simulations. 

Negative transmission times is a wave phenomenon, 
which together with superluminal propagation has been 
observed for light propagation through wave guides and 
through dispersive atomic media. It is not suprising that 
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FIG. 5: Wavepacket scattering on a BEC. Time series 
are shown for two different incoming momenta: At mean 
wavenumber koa = 9.75, the initial packet propagates to- 
wards the potential well containing the BEC and reflected and 
transmitted wavepackets appear. In accordance with Fig. ^ 
transmission is approximately 62% at this ko- At koa = 11.2 
we are in a transmission window and the reflected wavepacket 
has small amplitude and a double-peaked envelope. The fig- 
ures show both |t/(a;)| 2 and — |V(a;)| 2 , the particle contribu- 
tion and (— ) the hole contribution to the total density of scat- 
tering atoms. Note that the ground state quantum depletion 
also contributes to the density of atoms out of the condensate 
mode. This contribution is located inside the well and is not 
plotted here. 



the Bogoliubov-De Gennes equations show similar effects, 
and, indeed, Ray Chiao et al Jll[ have suggested a many- 
body interference mechanism that could lead to atomic 
transmission through condensates with negative trans- 
mission times. From our Bogoliubov analysis it is not 
clear if the result is due to this mechanism or if it is 
more closely related to the time delays, which may also 
be observed in normal wave packet tunneling through 
barriers Il2|. 




FIG. 6: Time spent in the condensate region. Curves are 
shown for both even and odd solutions. Where these curves 
approximately coincide, the retardation can be seen directly 
in wavepacket scattering. When they are very different, the 
translation to a time-dependent wavepacket formulation is 
less direct. Note the negative values of r e in the region around 
ka < 2: Wavepacket simulations show that here the peak of 
the transmitted wavepacket appears before the peak of the 
incident wavepacket has reached the condensate. For compar- 
ison, the two dashed curves show respectively the free motion 
and the sound wave tranmission times, 2a/v g and 2a/c. 



Note also that the analysis of both a real experiment 
and of wavepacket simulations is more complicated for 
massive particles than for light: As the transmission co- 
efficient depends strongly on k, there is a velocity filter 
effect, i.e., the transmitted wavepacket may move at a 
different speed than the incoming one. 



V. CONCLUSION 

It is well known that the Bogoliubov treatment leads 
to excited states ranging from phonon like disturbances 
of the mean field amplitude at low energies to particle 
like excitations at high enough energies. In the present 
case of a condensate trapped in a localized potential min- 
imum of finite depth, the eigenspectrum of low energy 
excitations corresponds to states which are particle like 
in some regions of space and phonon like in others. Like 
in normal scattering theory, wave packets are formed as 
superpositions, and they propagate as a consenqucncc 
of the phase evolution of the energy eigenstates. In our 
study, this propagation implies that particles are incident 
on the condensate, they propagate as phonons through 
the condensate, where they may be reflected back and 
forth between the condensate edges, and eventually they 
reemerge as reflected or transmitted particles. We have 
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determined the reflection and transmission probabilities 
and the phase shifts, which enable us to derive time de- 
pendent results from our stationary formulation. The 
main results of this analysis were double peaked distri- 
butions, reflected from the condensate and transmission 
with negative time delays. An experimental demonstra- 
tion of the latter phenomenon would be an interesting 
supplement to similar studies for light transmission. 

The particles are indistinguishable, and it is hence 
not meaningful to say that it is the same particle that 
emerges after the scattering process. The process, how- 
ever, is coherent, and intcrfcrometric experiments should 
be able to show that coherence is maintained, just as 
experiments with quantum correlated atoms should re- 
veal that also entanglement is faithfully preserved by the 
intermediate phonon excitation, in analogy with recent 
experiments where surface plasmons propagate quantum 



correlated photon pairs through sub wavelenght hole ar- 
rays jl3). We imagine that scattering experiments of 
the kind analyzed in this paper may be an ingredient 
in the study of controlled atomic dynamics, in particular 
for atoms and condensates trapped on chip architectures, 
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